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Otic o£ thfi ErOst widely used routinea in an algebraic min ip u La t ion 
system ia a polyDomial martipulfitien packageThe crucial 
operation in such routlrea ia chh eKtraction of the Greatest Cninmon 
IfivlBor (G^h) of twQ polynomials. This Operation is oriiclal because 
Of its fre^icric use and because it ia an espenaiT?® operation in regard 
to time aad space^ 

Experiments by Colllns^^^ have shown cJiat given two polynoTniale 
chosen at raodoiTip the GCH haa a high probability of being 1. taking 
into account this probEbility and the cost of cbtainiflg a GCD (some 
GCDs of polynoT.iials of degree 5 in tm or three vaii?ibles can take on 
the order of minute on the 70^4Ir appears that a quick method 
cf dcterminifig whether the GCO Is exactly ] would be profitable, Khlle 
no such coapLcte method ia known to exist, fi fail-safe procedure has 
been found nud la described here. A fati-safe procedure is cE^dtacteriz^d 
by the fseC that when It comes to a decision (in this caae that the GCD is 
1), then the decision is corrtet. However, the cooolusion (i,e, that tht 
CCD is 1) may be tt^ie, and tht procedure need not nrrivc at a decision 
regarding it. It is believed that the fail-aafe procedure presented; here 
(and its extension to the linear case) will arrive, at a decision quite 
frequently when the GCD ia actually 1, 

We shall first consider the case of polynomials of only one voriable. 
The extension to several vetiabLes will be made later, 

i.et P{x> “ x*’ 4 a^_j^ x" ^4 , , , +■ a^, integers, n^O 

Q(X] - + b^..^ k""-^ 4 4 by, b^ integers, m>D 

Let tlie operation of taking the GCD ef A,E (both fot integers and 
polynomials) be represented by CA,3]l , 



LflC {F{m) , d(K!l atld l*c 

P <x) - d (k) r (x)J Q Ck) = JCk) 3 (k>. 

fflr anj' ititejer a. then 

CPU). Q (a) ) • dU> 3 U> ) 

Shai; Cr(a) j s (a) ) n^ed noC te 1 aLchou")i 
(t (i) , a (x) ) ia 1 wtien the GCD operation is over poIynorLieLa ^ Note fuifther 
that the ahaalute val'je of d (a) is alvays Less than or e^^ual to the OCD i}f 
PCa) and <J(a) . 

The plan for the procedure is as foLLows: We would lita to substitute 
a fiumher. a. gay. into both PCx) and ()Cx) and ^et the GCD of the resulta. 
This integer a is se chosen that the leading term of dCx) will dominate the 
test of the polyaomtal, Ln faot, wo would Like the leading term of dfic) to 
he at least bvioc as LatgO as Che tost of dC'i<)r Thus i£ the rest of 

the terms of dC;:) are of opposite sign froc che leading termn we will be 

able to guarantee that che absolute value of dU) is at least one-half the 
absolute value of its Leading term. If dCK) be of degree k, then 

d(a) would be Least If (P(a>, Q(a) ) ia leaa than ^ja^* then d(x') 

cannot be of degree k. If (PCe). QCs) ) 5* Jja . then we do not have a 

decision regarding k. 

In order tO find a SuitAble candidate for a^ we will need some 
informat Lon about d.(x) . Without CK tensive calculation j tc appears that 
WE do not know the degree of d(j:) and certainly not its coefficients. We 
could, however, fiod a bound on the zeros of d(:t) since such zeros axe roots 
of both P(k) and 

If we let R be the rninlmutti of ft end K , where R and R axe bounds for the 

p q p q 

roDta of Pfs) and QCx) rospeceivtly, ehen we find the fallowing theorem 
applicable I 

Theorem Let RC^l) bound the absolute value of the roots of where 

dCx) is a polynomial of degree k>0. Thto for KS’lkR the absolute value of 
the leading term of dCx) will be grUAter than twice the absolute value of 
the remainder of d(K)r 
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Thus If we choose a>3feS , whexE k is the cilnlatLiffl flf m af,d n, 
theci w* could defjide that the CCD is a eonsUflC Cl.e. is af dt^rsE 

0) if (PCaJ, Q<aO } < w* ecmld decide that the OCD ia 1 if we 

prev-LQiisly had taken care of oonstatit GCDs fcy dlviains P[jc} atld 
ly the CCD of their coeffioientsJ 
Proof 

Given th&t the roots are boiinded hy R. the oondltlon regarding the 
leidLog term is harlest to satisfy when the coefficients of the reraainder 
Of the polyncffliftL are as Large as possible and of equal sigo. This OtOurs 
in tbe podynomtal , The oondlttcrl regarding the leading tera Gar. 

now he stated as 

tK+P)^ - ), x>0, R>0 

U+^]^ 



X + R 





1 


1 


R 


The fOllQ-vTing lemma will eoimplete the proof r 
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k a positive integer 
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By the foraula for flalte SMniS cf geometric progresslona 
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We have yet to ottaiii a cao^lete proc-ed'jre aince vfe laclj a bound 

on thfl roots or P{k) and Q{x}, A sinjlt bound for bh* roots of p(kl) 

is n Qiajc [ ^ 0.^ i4 n ■- 1} 
nn 

Tbls i* a onidt toiuid since tor (k+K? it yields JiR” instead of E. 

Several Other Oasil-y oalOLLLable bounds art hno^rn, but ve shall leavE this 
pxohleui to nunerLoal analysts. 

Now let us sun up the pLr-ocedttre defirue-d abovo. Wc find tho mininuiL 

of the bound on the absolute value of tbe roots of P(jc) and Q{jt). Call 

it* H* say. We pick an integer* a, say, such that a where k is the 

nlninuzD of the degrees of pfxj and Q.(x) . If we let g = (P|a}» Q[&)) then 
If g-^ ^ and if we previously ascertitned tb&t tho ooeffioients of P(x) 
and had nO' ectCKOTi OCD* then ve can guarantee that d{s) = {F{a) )= 1, 

Otherwise no decision is Hidft. 

Extensions to the Easic FroCedurt . 

Llntar PCD 

Let g * [PCa) . li(fl) ) (wh ere a is ohoaen by the procedure above). If 

S< ^ , ’h'e know d(k) ■ In Suppose, however* ^ f g < ^ at worst dCxJ COUld 

be linear in 3i, tbit ts, ■ CX + d (cft O). We propose now to extend the 

scheme defined above in or^er to detenaine utiether d(x) ia linear. If d(x} 
is linear then this extended procedure will find It; if it is not, then d(x) 
must bn 1; AQd then we have eKt&hdi?d th& range of applicability cf the 
proeedufs for determinitig whether the ■5CD la 1. 

As it turns out we can attempt to find all possible linear G':]D hy 
an extensive search since there are restrictive condtuians on o and d- 
We know that 

c|(a^, b^), dl {fc^i t^) where e|f 


means that s divides f 





w* be able to avoid tegtiti^ SOiiie er all possible Linear 

by maliiii^ use of tbe fsct that dCa) dividaa g. Fot any palrj 

Cjd, sneb that ac d does aot divide g, ue tjlti tondude that cm + d 

1,S net a possible d(K). If at l-d does divide g, then ue shall evaluate 
C c 

<}(- Both of these Tluj;bers ctust be zero. Tf they are^ 

d(xi = CK + dn If no possible pate of values for C and d works^ then 
dCx) = 1. 

Extension to Severi^l Variables 

One can estend the procedure for several variables cerely by 

substituting small integers for all but oc^e of the variable and 

deternining whether the resulting poJync-inials have a GCD of J. If 

they do the original polynomials did also I if they do not th&h ttP 

decision is ciade. Tills Idea is due to >9, >janove who progratTicd it in 
3 

MATHhAti ^ end existed prior to our result, 11 is likely that Mnnove'g 
idea by itself will yield a very large saving of effort since the 
tiDm|iutAtion of CGhs grows V*ty badly when the nnabar of variables 1? 
inoreAsed, 

The procodure was attempted on the pair of q^uartics given in ^l] , 
The evaluation was atteripted for values of a bet^’oen IGO and ISOr Host 
GCDs resulted in nucibers lass than thus yielding the conclusion of 
a GCD of 1, In Che ten ceigea that did not, tho linear Case settled the 
matter quickly since ec + d, for possible vEiloOs of c and d, did not 
divide the GCD with one exception at a ^ lO?', 

For J.arger pairs of polynomials (say, of degree 30) the evaluation 
will have Co be attempted on relatively large numbers' ton the ord&r of 
1000), Thus one would have to regort to variable precision arithmetic- 
We must consider, however, that the evaluation of the GCD o£ S^icb 
polynomials by strsiglit forwai'd methods wpuld also require variable 
precision integer o'rttbE.'ietlc. 
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It is hflTd to eBtiTEate> at pCBSent. Ch& effect oti running Clme 
that t)ie Intm^Liction af a fall-Eaj:E aethod vill tiai^re the CSLcuIaClan 
of a GCDe. a gain In ifeeil of an order of magnitude or better is not out 
of the flueation. If the fail-Eaft tethrtiilue gains aoceptanoe, then research . 
ii llhely to yield various metlioda with differing cosc# and ranfieg of 
applicability+ 
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Ev&lLi&l;ica , 

The guccega of the method de^sadsfn tfl Large ceaaure, on the likelihood 
that the sf r(a) and i(,a) will he fairly low. Below we give aa tnto’jr- 

agtnj result which shciws that if rfa) and E{aJ were chosen at random^ using 
a flat •dlstrihutioci,. then the probability that thelT CCD 1 e greater than 1 
li less than 0,4, the raeaningfulnese of this result in practical situations 
is hard tn gauge at preaent. It does, hotseyei:^ give nna great hope for the 
au-ccess of the method. 


Theorem , Let b he S-hy integer chosen aC random from 1,2 3 .h . 

The prohahility that gn integer less than b is relatively prime to b 


approacheg ^ 




as N appteacher 


Proof ; 

It la 


The number of integers less than b and relatively prime to 


b H (1 - —) j where p|b and p is a prime 


Tliat number can be written as 


b H (1 - Pr where PrCp|b) is the probability 

that p divides b which is 1 £f p]b and 0 otherwise, 

1 

tfow If we allow b to varyn Pr[p|b) = — , that is, the probability that b 

p + 

is even Is ^ * that b is divisible by 3 ia -j ^ etc. 

Thus, the nuEcber of integers lees Chen b and relatively prime to it 
becomeg * 
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b R fl -- bH (1 - where p ±s a prime. 




But^ 


^ U ’ fa) - 1 


Pt ^ Priffle n=l 

and 


h-1 


a = 


TT' 

% 


Thus 


bit (1 - ^2) - t I 2 

(The thesireiii used above may be found iu any bflob on the Bletann Seta TunctlOftJ , 

This astiraate for the probability is already known. The proof above la 
due Co W, Henneoan. 

Wa would Itally be Interested ih knowing the probability thiC 
Cr(a)j sfs)) t e, far email InC'e^et vsluaa of t, Certain theoretical 
ConslderaClorts indicate that for e = 10, Che probability is greater than 
0.9. However, Cheae considerations are based on cartain assumptions re¬ 
garding the distribution ef numbers relatively pTiz::e to a given number. 

In eKperiraanta we performed these estimates were off by one to two per cant. 
Nonetheless all signs point to the fact that if rfa) and s[a) were chosen 
at random, a high VSlwe for their GCD is very unlikely. 
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